Rabi oscillations of matter wave solitons in optical lattices 
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Inter-band Rabi oscillations of gap soliton matter waves induced by time dependent periodic forces 
in combined linear and nonlinear optical lattices are for the first time demonstrated. It is shown 
that under suitable conditions these oscillations can become long-lived. By switching off the external 
force at proper time it is possible to create either pure (stationary macroscopically populated) gap 
soliton states or linear combination of two gap solitons with appreciably long life-time. 
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I. INTRODUCTION 

In spite of its long history, the phenomenon of quan- 
tum tunneling still attracts a great deal of attention, 
particularly in connection with its macro-world mani- 
festations. In condensed matter physics tunneling phe- 
nomena occur in connection with inter-band transitions 
of electrons (from the valence to the conduction band) 
under the action of external perturbations such as elec- 
tric fields (Landau-Zener tunneling ^ ) or radiation fields 
(photons) of proper energy (inter-band Rabi oscillations). 
The phenomenon is highly facilitated by time periodic 
perturbations with frequencies matching the energy gaps 
between the bands. In the context of cold atoms Rabi 
oscillations (RO) between Bloch bands have been consid- 
ered both theoretically ^ and experimentally Pos- 
sible extensions of these studies to ultracold atoms and 
to Bose-Einstein condensates (BEC) give an interesting 
perspective for inter-band transitions and ROs detection 
at the macroscopic level. 

Macroscopic quantum tunneling phenomena are 
presently investigated in many physical systems, includ- 
ing superfluidity and superconductivity and, more re- 
cently, in optical systems with ultra-cold atoms Q and 
Bose-Einstein condensates (BECs) embedded in optical 
lattices (OLs) ((|, where the Landau- Zener tunneling has 
been recently observed f^. In contrast with solids, how- 
ever, BEC systems are intrinsically nonlinear, being de- 
scribed by the mean field periodic Gross-Pitaevskii equa- 
tion with the cubic nonlinearity arising from the inter- 
atomic two-body interactions. The correspondence with 
solid state physics occurs either at very low densities 
or when the interatomic interaction are "artificially" de- 
tuned to zero by means of external magnetic fields, using 
Feshbach resonances (see e.g. 

The presence of nonlinearity in the periodic Gross- 
Pitaevskii equation introduces modulational (dynamical) 
instabilities of Bloch wavefunctions 9] with substantial 
effects on the inter-band tunneling [7, 10]. The appear- 
ance of spatially localized states inside gaps of the un- 
derlying linear band is associated to these modulational 



instabilities. These states, also known as gap -solitons 
(GSs) first discovered in nonlinear optics [U\ and ob- 
served in fibers with Bragg gratings [l^l, in BECs were 
experimentally realized for the first time in p^ . Such 
states have energies which can scan the whole gap (or a 
large portion of it) giving rise to families of modes with 
different symmetry properties with respect to the OL (see 
e.g. [13 )• The existence of GSs with chemical potentials 
arbitrarily close to band edges introduces new possibili- 
ties for a stimulated Landau-Zener tunneling and RO of 
matter waves across the gap. 

The aim of the present paper is to show for the first 
time how the inter-band Rabi oscillations of matter wave 
solitons can be induced by means of combined linear and 
nonlinear OLs under the action of a time periodic linear 
potential. By switching off the external force at proper 
times we show that it is possible to put the system either 
in a stationary GS state or in a time dependent linear 
combination of two GS states which, in spite of the non- 
linearity of the system, can persist in this superposition 
for a reasonable long time. 

The paper is organized as follows. In section II we 
introduce the model equations while in section III we 
provide an analytical description of the long-lived Rabi 
oscillations of GS in terms of a modified two-mode model 
valid for GS states close to gap edges. A small detuning 
of the force frequency from the gap- width will be used as 
parameter to achieve the absence of matter leakage (emis- 
sion of radiation) from the soliton states. We show that 
under this condition the ROs become long lived and man- 
ageable. In Section IV we compare our analytical results 
with direct numerical simulations of the Gross-Pitaevskii 
equation. In this section we also discuss the possibility to 
create pure (stationary macroscopically populated) gap 
soliton states or linear combination of two of them which 
survive for appreciably long time, by switching off the 
external force sustaining the Rabi oscillation at a proper 
time. In the last section we design parameters for a pos- 
sible experimental observation of the nonlinear ROs and 
briefly summarize the main results of the paper. 
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FIG. 1: (a) Curves G±(V') separating domains where 
X±{y,G) > and XiCKG) < (above and below the re- 
spective hnes) with + (— ) denoting quantities evaluated at 
the top (bottom) edge of first gap. (b) Soliton norm A*' vs 
energy E of lowest symmetric (left curve) and antisymmetric 
(right curve) GS in the gap £- < E < £+ (the edges of the 
gap coincide with the boundaries of the i?-axis) for the pa- 
rameters V = 3, G = —0.5. Lower panels show shapes of the 
localized modes at points A - D of the panel b. 



II. MODEL EQUATION 

Although in the following we refer specifically to a 
BEC, the obtained results will be valid for all physi- 
cal systems modeled by the following periodic nonlinear 
Schrodinger (NLS) equation with time dependent linear 
potential: 

= -^p:cx+K{x)^ + J'{t) xi; + g{xM^^P, (1) 

Here l{{x) and g{x) are 7r-periodic linear and nonlinear 
lattices, respectively, and J-{t) is the time periodic ampli- 
tude of a linear potential. In BEC experiments the non- 
linear lattice can be created by optically [l^ or magnet- 
ically [l^ induced Feshbach resonances (for the sake of 
brevity we speak about OL) while the external force can 
be induced simply by a time periodic acceleration of the 
linear OL. In absence of the nonlinearity, g{x) = and 
of the Hnear potential J^{t) = 0, Eq. ([T]) reduces to the 
familiar band theory framework, with the Schrodinger 
operator 7i = —cP/dx'^ + 'U{x) and associated stationary 
eigenvalue problem TC(pnq{x) = £„{q)^nq{x), this pro- 
viding energy bands £„((/) and Bloch functions ipnq{x) 
with quasi-momentum q and band index n, normalized 
a-s \ifnq{x)\'^ dx = 1. We will concentrate on the low- 
est gap extending from the top of the first band, denoted 
as to the bottom of the second band, denoted as f+. 
The respective Bloch functions will be denoted by ip±{x). 

The condition for the existence of GS states at both 
edges of a gap can be expressed in terms of effective 
masses M± = {d^£±/ dq^^ and effective nonlinearities 
X± = r_^Q{x)\ip±{x)\''dx, as: M±x± < [13. In anal- 
ogy to what done in Refs. [l^jlillj we satisfy this condi- 
tion by considering OLs of the form: U{x) = — Fcos(2a;), 
g{x) = G + cos{2x), with V and G related to the physical 
amplitude of the linear OL V and to the scattering length 
as{x) = aso{G + cos{2x)) by V = VEn and G = {as)/asQ 



where {as) = tt~^ as{x)dx denotes the average of the 
scattering length over one period of the nonlinear OL. 
Here E^j = h^ir^ / {2md^) is the recoil energy, d is the lat- 
tice constant, a^o > is the amplitude of the scattering 
length modulation, and m the mass of an atom. Space 
and time variables in our formulas are measured in units 
d/TT, and h/Ej^, respectively, with the norm N related 
to the physical number of atoms TV by iV = '^°°!j^ A/" and 
with a± the transverse oscillator length. From Fig. [1^ 
we see that the OLs appropriate for our task are the 
ones with G,V taken in the region between curves G±. 
Notice that for fixed G and V a family of solutions, pa- 
rameterized by the soliton norm N = jV'prfa^ (for a 
cigar-shaped configuration bifurcates from the opposite 
band edges (panel b). The bottom panels of Fig. [1] il- 
lustrate change of the symmetry of a GS as the different 
edges of the gap. 

To stimulate RO between GSs we use a force T = 
v cos^ujt) with amplitude v 1 and frequency lo = Eg — 
2A slightly detuned from the gap width £g ^ £+ ~ 
(i.e. A <C £g). The oscillation must involve a nontrivial 
rearrangement of the matter at each cycle, due to the 
change of symmetry of the state and the corresponding 
dynamics can be sustained for a long times only if the 
force is properly designed. In the following we use A 
as a free parameter to find optimal inter-band tunneling 
conditions. 



III. TWO MODE MODEL FOR NONLINEAR 
RABI OSCILLATIONS 

An analytical description of the phenomenon can be 
made in term s of a two-level model. To this end we con- 
sider \fvju} ^ 1 as a small parameter and recall the 
semi-classical equation of motion dq/dt — — i/cos(wt), 
leading to the change of the phase of the wavefunction 



as follows e'[«(*)-'?(°5 



and take into 



account that in the leading order the solution is a super- 
position of the states bordering the gap edges: 



(2) 



Here A± = A±{x^t) ~ \fvT^ are slow varying ampli- 
tudes {dA±ldt ~ 5M±/5a;2 - {vl^)A±) of the Bloch 
states Lp± bordering the edges with energies /i± — £±=fA, 
providing [i^ — /i_ — lo and satisfying the above criteria 
A ~ J/ <^ w being small detuning towards the gap. Drop- 
ping details (see e.g. [TSIIl), we give the final system 
for the evolution of the amplitudes 



dAi 



1 

2M+ dx^ 



± AA± + 7Azp 



(x±l^±l'+2x|A^lVd 



(3) 



where x = fl^gix)iplf'idx and 7 = z; I-tv ^f-'^^- 
In absence of nonlinearity and for spatially homoge- 
neous amplitudes this system reduces to the model of 
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FIG. 2: (a,b) E versus time for v = (solid lines) and v = 
0.2 (dashed lines); (c,d) dynamics of the populations, pi,2(t), 
of the two bands for v — 0.1. Parameters are: V — 3, G — 
—0.644 and the initial conditions correspond to = —0.73 
(panels a,c) and E = 2.152 (panels b,d). Lower panels show 
the solitonic shapes corresponding to points A-F. In panels 
B and E soliton profiles are compared to stationary solutions 
(dashed lines) with the same A'^ of the initial conditions (A^ — 
0.61), but with energies in the vicinity of the opposite gap 
edge {E = 2.152 in panel B and £ = -0.73 in panel E). 



RO of two- level atoms in a light field [l9|. This con- 
firms the interpretation of the phenomenon as Rabi os- 
cillations of a gap soliton between the two energy lev- 
els fi±. Neglecting the second derivative and nonlinear- 
ity in Eq.Q, the inter-band tunneling time can be es- 
timated as half of the Rabi period, i.e. Tt„„ = ti'/^^, 
with fl the usual linear Rabi frequency fl = 1\p^ 



'7- + A2. 

By approximating Bloch function at the first gap edges 
as '/5+(a;) « tt"-^/-^ sin(x), ^-[x) Tr"-^/^ cos(a;), we 
readily estimate Tt„„ 



7r/(2^jy2/w2^_^2)^ The two 
GS states at gap edges can be estimated from the soli- 



ton solutions of Eq. jS]) given by A-^' ^ 0, A^'^ = 

with M±17± < 



XT 

and 17 



-»(0:FTA)tsg(,h[^-2A/zpl7^a;], 
fixed by N 



zp ii^cu ly - ^2|Mzpr2=p|/7r|AfzpXTl- In the 
following we use this solution to find an optimal design 
for long-lived ROs. To this regard, we remark that long- 
lived ROs can be achieved only if a minimal loss of mat- 
ter (ideally zero) from the solitons occurs during each 
oscillation cycle. For this one must impose the conser- 
vation of the energy and the conservation of N (com- 
plete exchange of particles between components at the 
turning points). Moreover, since the states must change 
symmetry every half cycles (being states at opposite gap 
edges) one must also require that the spatial localization 
of the two GSs at the turning points is approximately 
the same (this facilitates the symmetry exchange). The 
last condition is satisfied if M+x-i- = ^-X-^ while the 
conservation of iV and the conservation of the energy 
imply M_r2_ = M_(-i7+ and A = (ri_ — i7+) /6, respec- 
tively. Since f7± depends on N ^ M± ,x±, one can adjust 




FIG. 3: Long-lived (a) and fast decaying (b) Rabi oscillations 
in momentum space |c2(g)p. (c) Maximal population of the 
second band as obtained from NLS numerical integrations 
in the interval uj £ [£g — {Q- — Q+),£g] around the optimal 
value ujo (dashed line). Other parameters are v = 0.1, = 3, 
G — —0.644 (panels a,c], and G — 0.4 (panel b). The energy 
of the initial state is i5 = 0.73. 



these parameters to have all the above relations satis- 
fied with the driving frequency loq — £g — 2A chosen as 
ujQ =fg — (0_— ri_(_)/3 for an optimal design (see FigjSJ;) . 
Notice that for u! ujq, while N will be conserved, the 
soliton energy will not be perfectly conserved during the 
tunneling process and that for fixed M± ,x±, a change 
of N requires the adjustment of the optimal frequency. 
Contrary to the linear case, where the optimal transfer 
between levels is achieved for zero detuning, here the op- 
timal transfer is obtained for A 7^ 0. This is a direct 
consequence of the nonlinearity of the system which re- 
places linear atomic wavefunctions with energies exactly 
at gap edges with nonspreading GSs wavepackets with 
energies detuned from band edges inside the gap. In the 
limit of zero nonlinearity {N 0) GSs tend to extended 
Bloch states at gap edges and the linear optimal condi- 
tion (A = 0) for RO is recovered. 



IV. NUMERICAL RESULTS 

In order to check these predictions we have performed 
direct numerical integrations of Eq.([T]) with initial con- 
dition in a form of a stationary GS close to the bottom 
edge of the gap. The frequency detuning, the parameter 
values of the OLs, and the initial conditions were chosen 
to met the requirements of conservation of energy, soli- 
ton norm and widths in the solitonic states during the 
oscillation, as described above. 



A. Nonlinear Rabi scillations and tunneling time 

The results are depicted in Figl2] where, in panels a,b, 
we show that minima and maxima of the periodically 
varying energy of the system well match the GS levels, 
what represents the most clear evidence of the inter-band 
tunneling. Notice that the frequency of the energy os- 
cillations increases with the force strength v and corre- 
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FIG. 4: (a) Tt^n fs v (in the logarithmic scales); (b) evolution 
of the soliton norm in the lowest band A'^_ sX v = 0.01, ob- 
tained from the direct numerical simulation of ([1} (solid lines) 
and from the two-mode approximation ([3]) (dashed lines). 



lates with similar oscillations of the band populations 
(Fig. [21;, d) computed as pi_2 = \c\,i{(l)'\^dqlN , where 
c-nisi) = J^^ip(j)n,qdx. One also observes change of the 
wave-packet symmetry: initially symmetric (with respect 
to a; = 0) soliton excited in the vicinity the lower gap edge 
(Fig. [2] A), in the vicinity of the upper gap edge trans- 
forms into an odd solution (Fig. [2]B). Similarly, initial 
odd soliton (Fig. [2] D) is transformed into the even one 
(Fig.[2]E). After one cycle of oscillations (points C and F 
in the upper panels of Fig. [2|) , the soliton approximately 
restores its initial shape (c.f. panels A and C, and panels 
D and F). 

To compare RO under optimal and arbitrary condi- 
tions, we show in Fig. [3K,b, the time evolution of the 
GS states projected on the second band (i.e. the popula- 
tion \c2{q)\'^) as obtained from numerical integrations of 
the NLS in the case of an optimally chosen nonlinearity 
G (Fig. [3^) and for a non optimal value of G for which 
small-amplitude GSs do not exist close to the upper band 
edge (Fig.[3jD). We see that while in the first case the ROs 
are long-lived, in the latter case the soliton is destroyed 
after few periods of oscillation. From Fig. [SJ; we also see 
that the numerical frequency for the maximal transfer of 
atoms from one band to another is in very good agree- 
ment with the analytical value loq estimated from the two 
mode model. 

In Fig. m we present the comparison of the tunneling 
time Ttun obtained from the direct numerical simulations 
and in the framework of the two-mode model. One ob- 
serves that Ttun decreases with the increasing the force 
strength i/, and increases with the amplitude of the linear 
OL. The discrepancies between the exact and approxi- 
mated models are negligible for moderate amplitudes of 
the linear OL (namely V <2 in Fig. UK), and become ap- 
preciable for a deep OL and small external force strength 
(F = 3 and 1/ < 0.01 in Fig. [4}i and = 3 in Fig. Hd). 



B. Superposition of two gap soliton states 

It is worth to note that by switching off the linear po- 
tential at a time tg = (2n -I- l)Ttun {n = 0,1,...) the 
state of the system will coincide with one of the station- 



FIG. 5: (a) The soliton energy E{t) for E{0) = -0.73 (solid 
lines) and E{0) — 2.152 (dashed lines), with the linear po- 
tential switched off at Ttun ~ 42.6 (the vertical dashed line). 
Parameters of OLs are fixed as in Fig. [2] but with i/ — 0.1. 
Panel A shows the soliton shapes at point A of panel a (solid 
lines), compared with the stationary solutions (dashed lines) 
depicted in Fig.[2j3. (b) Population densities |c±(t)|^ after the 
switching off of the linear potential ai t = 20.6. 



ary states close to a band edge (see Fig. [S^ and panel A 
in the same figure). On the other hand, from Eq. ([2]) 
we see that by switching off the linear potential at any 
other time tg ^ {2n + l)Ttun will produce a state which 
is a linear superposition of the two stationary states at 
the edges of the gap: tps = c^{t)^/j^ + c+{t)'ip+, where 
c± = J ip{x,t)ip±{x)dx are complex amplitudes and 
ip± are the numerically exact GSs with the same N used 
in direct simulations. It is remarkable that, in spite of 
the nonlinearity of the system, the decay rate of the su- 
perposition state (slopes of c±) may be small enough for 
the state to be observed. The decay rate depends, how- 
ever, on the soliton norm, and we find that for A'' = 0.17 
is ~ 4.5 • 10-^ at TV = 0.61 is - 3.5 • lO"'' (see Fig. ^) 
and for N — 1.24 it becomes ~ 1.6 • 10^'^, as consequence 
of the increased effective nonlinearity in the system. 



V. DISCUSSION AND CONCLUSIONS 

We discuss a possible parameter design for experimen- 
tal observation of long-lived Rabi oscillations between 
gap soliton states at opposite edges of a band gap. Re- 
ferring to a Li condensate in a trap with a± = 2 /im, 
d = 1 /xm and with (as) = — 1.288 nm, Oso = 2nm cre- 
ated by an optically induced Feshbach resonance, the RO 
depicted in Figs. [2t,d can be achieved by an external 
force of 1.45 • 10"^'' N (corresponding to an acceleration 
of the OLs « 124 m/s ), oscillating with the frequency 
20.38 kHz and applied to a GS with about TV = 980 
atoms (this corresponding to A^ = 0.61). The tunnel- 
ing time and the decay rate of the superposition state in 
Fig. [5t are estimated as 0.965 ms and ^ 15.5 s~^, respec- 
tively. 

In summary, we have shown for the first time how 
to achieve long-lived inter-band Rabi oscillations of gap- 
solitons matter waves in combined linear and nonlinear 
optical lattices. We also showed that by switching off 
the external force at proper time it is possible to create 
either pure (stationary macroscopically populated) gap 
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soliton states or linear combination of two gap solitons 
with appreciably long life-time. 

We finally remark that similar phenomena are ex- 
pected to occur also in arrays of optical waveguides pe- 
riodically modulated along the propagation direction. 
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